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Summary 


The heat transfer coefficients of air-glass, argon-glass, and 
argon-aluminum suspensions were measured in horizontal and verti- 
cal 0.870 inch, inside diameter, tubes- Two grades of glass beads, 
averaging 21.6 microns and 36.0 microns diameter by weight, were 
suspended in the gas at gas Reynolds numbers between 11,000 and 
21,000, at loading ratios between 0 and 2-5. The aluminum powder, 
averaging 30.0 microns diameter by weight, was suspended in the 
argon at argon Reynolds numbers between 11,000 and 21,000, at 
loading ratios between 0 and 1.0. The suspension was continuously 
recirculated in a closed loop system, originally designed for flow 
measurements but modified to permit heat transfer experiments - 

The presence of particles generally reduced the heat transfer 
coefficients below those of the clean gas at corresponding gas 
flow rates. Reductions of up to 44% and 52% were achieved in the 
vertical test section for the 21.6 micron and 36.0 micron glass 
beads respectively. Reductions of up to 38% and 44% were also 
achieved in the horizontal test section for the 21.6 micron and 
36.0 micron glass beads respectively. The exceptions to the 
general occurrence of Nusselt number reduction appeared either at 
low loading ratios or at low Reynolds numbers. The Nusselt number 
increases were always less than 20% and were usually less than 10%. 
With the exception of cases where deposition complicated the 
results, the suspension Nusselt numbers were independent of the 
carrying gas (air or argon) and were independent of particle prop- 
erties (other than diameter) . 

The Nusselt number reduction in the vertical test section was 
considered as a coincident result of the drag reduction that had 
been observed earlier using the same experimental apparatus. The 
thickening of the viscous sublayer, the cause of the drag 



reduction, was also believed to be the primary cause of the 
Nusselt number reduction. Visual observations of non-uniform 
particle concentration and of deposition in the horizontal test 
section indicated that the particles here took a more active role 
in the heat transfer mechanism. Therefore the Nusselt number 
reduction in the horizontal test section could not be attributed 
solely to the distortion of the fluid structure by the presence 
of particles. Thickening of the sublayer was probably more 
important in the vertical test section; and, particle deposition 
was probably more important in the horizontal test section. 

The most unusual type of deposition was exhibited by the 
aluminum powder. Rock-hard deposits appeared on any protuberance 
into the flow channel. The aluminum powder welded itself to 
metallic objects by plastically deforming on collision with 
initial non-uniform deposits formed by other means. The other 
means of deposition, which were also appropriate to the glass 
beads, could have been caused by electrostatic attraction, gravity 
sedimentation, or thermophoresis. 

Introduction 

This report represents the final phase of a three phase 
investigation (ref. 1) into the existence, causes, and effects of 
drag and Nusselt number reduction. The first phase consisted of 
friction factor measurements which confirmed the existence of the 
gaseous drag reduction that was reported in an earlier study ' 
(refs. 2 - 3) . The first phase also indicated that the behavior 
of a gas-solids suspension compared to that of the clean gas was 
a_ strong function of test section orientation, a moderate function 
of particle diameter, and a weak function of gas Reynolds number. 
The second phase of the investigation consisted of gas velocity. 
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gas turbulence intensity, gas turbulent spectra, and particle 
velocity profile measurements designed to elucidate the reasons 
for the drag reduction. These measurements (taken in a vertical 
test section) revealed that gaseous drag reduction could be 
characterized by thickening of the viscous sublayer in the 
presence of particles. The third phase of the investigation, 
reported here, consisted of heat transfer measurements with the 
aim of providing an analogy between drag reduction and Nusselt 
number reduction - 

Previous investigators have reported decreased heat transfer 
when particles were added to a turbulent gaseous stream. A 
summary of earlier experimental evidence of Nusselt number 
reduction is shown in Table 1. Although experimental 
restrictions may have dictated the results, Nusselt number reduc- 
tion has generally been limited to small particle sizes, small 
pipe diameters, low loading ratios, and low turbulent Reynolds 
numbers . 

Only two of the investigators. Depew and Boothroyd, have 
measured both drag and Nusselt number reduction in the same test 
apparatus. However, neither^Depew nor Boothroyd varied particle 
diameter, particle type, or test section orientation. In this 
phase of the investigation, three different particle diameters, 
two different particle types, and two different test section 
orientations were considered. The particles were two grades of 
glass beads, averaging 21.6 microns and 36.0 microns diameter by 
weight, and one grade of aluminum powder, averaging 30.0 microns 
diameter by weight. The test sections were vertical and hori- 
zontal 0.870 inch, inside diameter, tubes. The suspensions were 
circulated at gas (air or argon) Reynolds numbers between 11,000 


3 



Table 1 


Experimental Evidence of Nusselt Number Reduction 


Reference 

Solids Diameter 

Pipe Diameter 

Loading 

Reynolds 


and Type 

and Orientation 

Ratio 

Number 

Farbar 
et al. 1957 

10-210U 

0 . 69" vertical 

0-0.20 

13,500 

(re'f. 7) 

silica-alumina 


0-0.40 

19,000 




0-0.70 

23,000- 

27,000 

Tien 

et al. 1962 

30u glass 

0.71" vertical 

0-2.5 

15,000 

(ref. 8) 



0-3.5 

30,000 


30|J lead 


0-3.0 

15,000 




0-3.5 

30,000 


200ia glass 


0-3.5 

15-, 000 




0-3.5 

30,000 


200|a lead 


0-3.5 

15,000 




0-3.5 

30,000 

Abel 

et al. 1963 

l-lOp graphite 

0.89"vertical 

0-5.0 

22,000- 

(ref. 9) 




65, 000 

Depew 

et al. 1963 

30u glass 

0.71" vertical 

0. 5-3.0 

13, 500 

(ref. 10) 



0. 5-3.0 

27,400 

Farbar 
et al. 1963 

30|j glass 

0.67" vertical 

0-1.5 

15, 300 

(ref. 11) 



0-1.7 

26,500 


70d glass 


2-3.0 

15,300 

Jepson 
et al. 1963 

211-295^ 

1 . 5" vertical 

0-2.0 

14, 300 

(ref. 12) 

sand 


0-4.0 

21,450 




0-6.0 

28,600 




0-11.0 

42,900 


76-104P 


0-5.1 

14, 300- 


sand 



28,600 


422-599P 


0-10.1 

28,600- 


sand 



71,500 
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Table 1 (cont.) 








Reference 

Solids 

and 

Diameter 

Type 

Pipe Diameter 
and Orientation 

Loading 

Ratio 

Reynolds 

Number 


Jepson 
et al. 1963 

850-1200P 

1 . 5" vertical 

0-2.5 

21,450- 

(ref. 12) 

sand 



71,500 

(cont. ) 

1200-1680|i 


0-8.5 

21,450- 


millet seed 



71,500 

Garret 1964 

243|J. sand 

3"vertical 

0-3.75 

23,000' 

(ref. 13) 

159p sand 


0-5.5 

30,700 

Szekely 
et al. 1966 

150, 300, 600, 

cyclone 

0-3.0 


,(ref. 14) 

1200|i sand 
200p bronze 
175p iron 




Sukomel 
et al. 1967 
(ref. 15) 

65p graphite 

0 . 30"horizontal 

0-8.0 

15,000 

Boothroyd 
et al. 1970 

0-40p zinc 

1" vertical 

0-20. 

35,000 

(ref. 16) 



0-8.0 

53,000 




0-8.0 

80,000 




0-5.0 

100,000 



2"vertical 

0-3.5 

35,000 




0-3.0 

53,000 




0-4.0 

80,000 




0-2.5 

100,000 



3" vertical 

0-2.0 

35,000 




0-2.0 

53,000 




0-3.5 

80,000 

Depew 

et al. 1970 

30|J glass 

0 . 71"horizontal 

0-4.5 

10,000 ■ 

(ref. 17) 



0-7.5 

15,000 


- 


0-4.0 

30,000 


200p glass 


0-6.0 

15,000 




0-5.0 

30,000 
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and 21,000 and -at loading ratios between 0 and 2.5. 

The range of parameters chosen here for the heat transfer 
work was identical to the range of parameters considered earlier 
during pressure drop work with a few exceptions. During pressure 
drop measurements, 15.0 and 55-0 micron diameter glass beads were 
also considered. However, the 15.0 micron diameter glass beads 
were found to be extremely difficult to circulate as their small 
size permitted rapid deposition and plating on the walls of the 
closed loop. The 55.0 micron diameter glass beads were found to 
degrade after prolonged passage through the pump used to circulate 
the suspension in a closed loop operation. The 15.0 and 55.0 
micron diameter glass beads would therefore not provide meaningful 
information after the long warmup periods required for heat trans- 
fer measurements. 

The aluminum powder, not used during pressure drop measure- 
ments, was selected for heat transfer tests because it had differ- 
ent thermal properties from the glass beads and might give addi- 
tional information on the heat transfer mechanism. Argon, also 
not used during pressure drop measurements, was chosen as the 
conveying medium for circulation of aluminum powder, because it 
would eliminate any possible hazard in the closed loop system by 
the buildup and release of static electric charges. 

Experimental Apparatus and Procedures 

The suspension was circulated through a closed loop system, 
originally designed for pressure drop and flow measurements but 
modified for heat transfer measurements. The closed loop system 
as modified for heat transfer is shown in Figure 1. The loop 
could have been described as a 12 foot high, 18^ foot long 
rectangle generally constructed of 1 inch outside diam.eter (O.D.), 
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0.065 inch wall, type 304 stainless steel tubing. The motive 
force for flow was provided by a centrifugal circulator which was 
capable of pumping gas and solids simultaneously. 

Two primary test sections were provided, a heated vertical 
section and a cooled horizontal section. The vertical test 
section was fabricated from a 30 inch length of 1 inch O.D., 0.065 
inch wall stainless steel tubing and located 92 tube diameters 
upstream of the nearest flow disturbance, the reinforced elbow 
near the circulator outlet. Miniature chrome 1 -a lumel thermo- 

r 

couples were located at either end of the test section and posi- 
tioned to measure centerline bulk temperature. The thermocouples 
were inserted through mounting bushings welded to the tube wall. 
The bushings seated securely on the 1/16 inch O.D. thermocouple 
sheath and did not permit leakage^ Two diametrally opposed 
chrome 1 -a lumel thermocouple junctions were located every six 
inches along the test section to measure outside wall temperatures. 

The heat input in the vertical test section was supplied by 
an electrical resistance tape wrapped around the tube. The tape 
was controlled by a 0 to 120 volt range voltage regulator and 
supplied a maximum of 384 watts. Several inches of fiberglass and 
asbestos wool insulation were wrapped around the tape to minimize 
heat conduction to the atmosphere. A final layer of aluminum foil 
was added to minimize radiative heat losses. 

A 48 inch long, 2 inch O.D., copper, water cooling jacket was 
located 78 inches downstream of the first reinforced elbow and 
served to identify the horizontal test section for heat transfer. 
Chrome 1-alumel bulk thermocouples were located at either end of 
the test section to measure suspension temperatures. ASTM A9F 
thermometers, reading to 1°F, were inserted through nylon fittings 
in the water jacket to measure the inlet and outlet water 
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temperatures. Another 48 inch long, 2 inch O.D., copper, water 
cooling jacket was located 72 inches from the unreinforced elbow 
in the lower horizontal section. Although this section was 
equipped with the same instrumentation as the upper horizontal 
section and was occasionally used as a test section, its main pur- 
pose was to protect the pump from excessive fluid temperatures. 
Five chrome 1 -a lumel thermocouples on the circulator were used to 
monitor housing and bearing temperatures. All thermocouple 

readings were taken from a 24 point dial recorder. The dial read 

O O' . o 

to 5 F over an 800 F range but could be estimated to 1 F . 

The gas flow rate was measured with a sharp-edged orifice 
even with particles in suspension. Earlier studies (refs. 1-3) 
had shown that the calibration of an orifice was unaffected by 
dilute volumetric concentrations of solids suspended in the gas. 
The solids flow rate was determined with the aid of a cantilevered 
target meter. The operation of the meter was based on the fact 
that the deflection of the cantilever (determined by strain 
gauges) was directly proportional to the force of the fluid 
impinging on the target. The force of the gas alone was known 
from the flow rate given by -the sharp-edged orifice. The force of 
the solids was related to the solids flow rate. During calibra- 
tion of the target-meter, the solids flow was determined directly 
in an open loop operation and used to separate the response of the 
target-meter to solids from the response of the target-meter to 
gas alone. Once the separate responses of the target-meter to the 
various particle sizes and to the gas was known, the target-meter 
was used to calculate the loading ratio in closed loop operation. 
The detailed calibration of this device will not be given here as 
it is lengthy and is described elsewhere (ref. 1). 

Solids were added to the loop before a set of data runs by 
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pouring a given weight of sifted particles into the solids inlet 
port located near the circulator intake. The port was sealed and 
the circulator was started to disperse the particles through the 
loop. Not less than a one-hour warmup period with the circulator 
running preceded each set of heat transfer runs with a given 
amount of particles in the loop. The cooling water temperature 
and the maximum wall temperature in the vertical test section had 
to stabilize before reliable data could be taken. At the begin- 
ning of each set of data runs, the cooling water flow rate (moni- 
tored on a 0 to 1 gallon-per-minute rotameter reading to 0.02 gpm), 
the circulator shaft speed, and the voltage regulator setting were 
set to the desired values. Because of the limitations on the 
allowable tape temperature, the voltage regulator setting was 
adjusted so that the maximum wall temperature in the vertical test 
section did not exceed 420°F. Although a one-hour warmup was 
required between sets of runs, only a fifteen minute period was 
required between individual runs where only the gas Reynolds num- 
ber and not the amount of solids in the loop was varied. For 

argon gas measurements, the system was purged of air by forcing 

/ 

the denser argon into the lower loop and venting the lighter air 
from the upper loop. Additional details of the experimental 
apparatus and procedure are described in reference 1.' 

Clean Gas Test Results 

Prior to suspension tests, heat transfer measurements were 
used to calculate the Nusselt numbers for the two gases, air and 
argon. The direction and instrumentation of the heat transfer was 
different for the vertical and horizontal test sections. Two 
different techniques were therefore used to calculate the gas side 
heat transfer coefficient and Nusselt number. 


9 


from: 


In the vertical test section, the heat input was calculated 


Q 


W C 

g pg 


(T -T 
OUT IN 


) 


( 1 ) 


The heat input was equated to: Q = UA(T^ - T^) (2) 

where T^ was an average wall temperature weighted in accordance 

with the linear spacing of the 12 thermocouples along the heated 
section with UA given by: 


UA 


'ss 1 

+ 


1 TV 

k A 11 
ss 


(3) 


With the heat input and all the temperatures as known quantities, 
the gas side heat transfer coefficient was calculated directly. 

In the horizontal test section, a counterflow heat exchanger 
equation was used to evaluate UA. The gas side heat loss was 
calculated from: 


Q 


w c 
g pg 



T ) 
OUT 


(4) 


The heat loss was equated to: 


(T -T ) - (T -T ) 
' OUT WOUT ' IN WIN 


Q = UALMTD where LMTD= 


Ln 


r^OUT_^WOUT^ 
vin"\in / 


and UA = 


ss 

k A 
ss 


+ 


h A 
o o 


h.A. 
1 1 


(5) 


In operation, the water temperature remained almost constant from 
the inlet to the outlet of the cooling jacket. The logarithmic 
temperature difference was almost exclusively dependent on varia- 
tions in the gas temperature. 'The water flow rate was high enough 
that UA was sensitive only to the gas side coefficient. The 
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resistance of the stainless steel tube wall was negligible in both 
test sections. 

The Nusselt number was evaluated at the average bulk gas 
temperature. The results for air and argon are shown in Figure 2 
where Nusselt number is plotted versus Reynolds number. The 
measured Nusselt numbers were about 10 percent lower than the 
Nusselt numbers calculated from the Sieder-Tate equation which 
correlates experimental data within wide 30% limits (ref. 4). 
Prandtl number dependence was not an important factor because the 
one-third power of the Prandtl number for air and argon is 
essentially constant for the experimental range of temperatures. 
The clean gas values were the same in both the vertical and hori- 
zontal test sections and the lower horizontal cooling section as 
well. Clean gas Nusselt numbers were measured before and after 
suspension tests and were found to be unchanged. 

Suspension Test Results and Discussion 

The suspension heat transfer, measurements were used to calcu- 
late Nusselt numbers for the two gases carrying suspended glass 
beads or aluminum powder. The circulation of suspensions required 
revision of the clean gas formulae for heat input in the vertical 
test section and heat loss in the horizontal test section. The 
mass flow rate now included the contribution of the suspended 
solids: 

W = W (1 + ri) (6) 

s g 

The specific heat of the suspension was calculated from a 
weighted average of the separate gas and solids specific heats: 
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( 7 ) 



The final calculations and correlations of Nusselt numbers were 
based on the gas properties (density, thermal conductivity, 
specific heats, and viscosity) in order to provide a direct 
comparison with the clean gas results. For the loading ratios of 
this study, the assumption that the effective thermal conducti- 
vity of the suspension was the same as the gas thermal conducti- 
vity introduced negligible error (ref. 5) . 

The results of the vertical test section heat transfer 
measurements are shown in Figure 3 for suspensions of the 21.6 
micron glass beads in air and in Figure 4 for suspensions of the 
36.0 micron glass beads in air and argon and suspensions of 30.0 
micron aluminum powder in argon. The graphs were plotted as the 
ratio of the suspension Nusselt number to the clean gas Nusselt 
number versus loading ratio at corresponding gas Reynolds numbers. 
Nusselt number reduction was indicated by values of the Nusselt 
number ratio below unity. Figure 4 shows the envelope of the 
actual data. The upper limit generally represented the argon- 
glass results. The lower limit generally represented the air- 
glass and argon aluminum results. 

The glass-air suspension measurements indicated that the 21.6 
and 36.0 micron glass beads were much more effective in reducing 
Nusselt numbers than in reducing friction factors (ref. 1). 
Reductions of up to 44% and 52% were achieved in the vertical test 
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section for the 21.6 and 36.0 micron glass beads respectively. 

Like the friction factor results, there was a minimum loading 
ratio for the onset of Nusselt number reduction. Also like the 
friction factor results, the onset of Nusselt number reduction was 
further delayed with increased particle size. Neither of the two 
glass bead suspensions showed any Reynolds number effect on the 
onset of Nusselt number reduction or on the amount of Nusselt 
number reduction. 

The suspensions of 21.6 micron glass beads showed a monotonic 
decrease in Nusselt number ratio after the onset of Nusselt number 
reduction. Before the onset of Nusselt number reduction at a 
loading ratio of 0.12, a moderate increase in Nusselt number was' 
noted. The increase was attributed to an initial widely dispersed 
deposit on the test section wall. The deposit was apparently 
insufficient to cover large heat transfer areas but was sufficient 
to disturb the viscous sublayer, producing additional turbulence 
and heat transfer. 

The suspensions of 36.0 micron glass beads in air showed a 
tendency toward a flat minimum at a loading ratio of 1.2. Before 
the onset of Nusselt number reduction at a loading ratio of 0.4, a 
slight increase in Nusselt number was noted. The increase was 
also attributed to a widely dispersed deposit on the vertical test 
section wall. That the increase was smaller for the 36.0 micron 
glass beads than for the 21.6 micron glass beads was consistent 
with the relative ability of the two bead sizes to plate out. For 
both bead sizes, the deposit, rather than accumulating further, 
was probably worn away by heavier loadings and did not affect the, 
vertical test section results past the onset of Nusselt number 
reduction. 

The measurements of Nusselt nuitbers for suspensions of the 
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36.0 micron glass beads in argon were performed in order to deter- 
mine the effect of different gas properties on the Nusselt number 
reduction. The argon appeared to delay the onset of Nusselt 
number reduction to a loading ratio of 0.6 compared to 0.4 for 
air. However, the maximum amount of Nusselt number reduction was 
approximately the same for both gases. The general shape of the 
Nusselt number ratio versus loading ratio curves was also the same 
for both gases. These results were not surprising since the 
Prandtl numbers of air and argon were similar in the working tem- 
perature range. 

The heat transfer measurements for suspensions of aluminum 
powder in argon indicated that particle thermal properties had no 
effect on Nusselt number reduction in the vertical test section. 
The 30.0 micron aluminum powder gave similar results as the 36.0 
micron glass beads. The similarity in behavior of the two differ- 
ent particles in the vertical test section suggested that the 
Nusselt number reduction was a coincidental result of the drag 
reduction mechanism; that is, a thickening of the viscous sub- 
layer. 

The results of the horizontal test section heat transfer 
measurements are shown in Figure 5 for suspensions of the 21.6 
micron glass beads in air and in Figure 6 for suspensions of the 
36.0 micron glass beads in air and argon and suspensions of 30.0 
micron aluminum powder in argon. Significant Nusselt number 
reduction, up to 38% and 44% for the 21.6 and 36.0 micron glass 
beads in air respectively, was also achieved for the horizontal 
test section. However, the horizontal test section Nusselt number 
results were not as consistent with the friction factor results 
as those in the vertical test section (ref. 1) . (In the vertical 
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test section Nusselt number reduction occurred under the same 
conditions as the previously observed drag reduction. In the 
horizontal test section, Nusselt number reduction occurred even 
for those cases where drag increases and not drag reduction were 
previously observed.) Gravity sedimentation and thermophoresis 
causing deposition on the cool wall of the horizontal heat trans- 
fer test section apparently complicated the results. 

The 21.6 micron glass beads (like the previous friction 
factor results) showed very similar results in both horizontal and 
vertical test sections with the exception of a detectable Reynolds 
number effect in the horizontal test section. The amount of 
Nusselt number reduction progressively decreased with decreasing 
Reynolds number. The onset of Nusselt number reduction was pro- 
gressively delayed, from loading ratios of 0.07 to 0.16, with 
decreasing Reynolds numbers until the lowest Reynolds numbers, 
about 12,000 were reached. At Reynolds numbers of about 12,000, 
Nusselt number reduction was not generally noted. The observed 
loss of entrainment at these Reynolds numbers could presumably 
explain the anomalous results as the presence of large quantities 
of solids in direct contact with the tube wall would drastically 
change the heat transfer mechanism. With the exception of the low 
Reynolds number results, the trends of the heat transfer measure- 
ments . for the 21.6 micron glass beads were reasonably consistent 
with the trends of the previous pressure drop measurements. 
However, it could not be concluded that the controlling mechanism 
in the horizontal test section was the same for both cases. 
Thickening of the sublayer was probably more important for the 
pressure drop measurements; and, particle deposition was probably 
more important for the heat transfer measurements. 

The results for suspensions of 36.0 micron glass beads in air 
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were less predictable and systematic than the results for the 
21.6 micron glass beads. Whereas, friction factor increases were 
previously noted at all test conditions for the 36.0 micron glass 
beads (ref. 1) , Nusselt number increases were found only for the 
loading ratios less than 0.4 or for Reynolds numbers of about 
12,000. In general, substantial Nusselt number decreases were 
observed. The friction factor and heat transfer results for the 
36.0 micron glass beads tended to confirm the conjecture that 
particle deposition was more important for the Nusselt number 
reduction in the horizontal test section than any uniform sublayer 
changes. The onset and amount of Nusselt number reduction did not 
follow the consistent pattern observed with tests of 21.6 micron 
glass bead suspensions. It was apparent that the Reynolds analogy 
which predicts corresponding effects on heat transfer and momentum 
transfer was not valid for horizontal flow. The particles did not 
just distort the structure of the turbulence as they did in verti- 
cal flow, but actually became an active part of the heat transfer 
mechanism. 

The data obtained with suspensions of 36.0 micron glass beads 
in argon were in agreement with the air suspension data. However, 
the aluminum-argon suspension data always indicated Nusselt number 
increases regardless of loading ratio or Reynolds number. These 
results were probably caused by deposition of the powder on the 
cool wall of the horizontal test section. The unusual nature of 
the aluminum deposits probably disturbed the flow pattern and 
produced additional turbulence. Reductions in local heat transfer 
rates directly at the deposit locations were probably quite small 
because of the high thermal conductivity of the aluminum. There- 
fore, it was likely that the additional turbulence/ aiding heat 
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transfer^ more than offset the local fouling factor Sy hindering heat 
transfer. 

The form of the aluminum deposit was quite unusual, espe- 
cially on protuberances into the flow channel. A rock-hard build- 
up of aluminum was observed on any surface subject to direct bom- 
bardment. Very heavy accumulations were noted on thermocouples, 
the target meter, and the internals of the circulator. Figure 7 
shows the type of deposit that was found on all thermocouples 
after bombardment by aluminum powder. 

The propensity of the aluminum powder to bond so solidly to 
the loop walls and equipment may have been the result of its wide 
size distribution. Although the aluminum powder averaged 30.0 
microns diameter by weight, the powder included particle sizes 
ranging from 5.8 to 55.0 microns diameter. Boothroyd (ref. 6) 
maintains that a polydisperse powder will afford a much greater 
number of additional cohesive contact points than a narrowly 
sized powder (like the glass beads) . The aluminum powder welded 
itself to metallic objects by plastically deforming on collision 
with initial non-uniform deposits formed by other means. The 
other means of deposition could have been caused by electrostatic 
attraction, gravity sedimentation, or thermophoresis. 

Conclusions 

The heat transfer coefficients of air-glass, argon-glass, and 
argon- aluminum suspensions were measured in horizontal and verti- 
cal tubes. The results of these tests led to the following con- 
clusions: 

1. The presence of dilute loadings of small particles in a 

turbulent stream reduced the heat transfer coefficient below 
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that of the clean gas at corresponding gas flow rates. 
Reductions of up to 44% and 52% were achieved in the vertical 
test section for 21.6 and 36.0 micron diameter glass beads 
respectively. Reductions of up to 38% and 44% were also 
achieved for the horizontal test section. 

2. The Nusselt number reduction in the vertical test section was 
attributed to a postulated drag reduction mechanism observed 
previously under the same experimental conditions; that is, 
thickening of the viscous sublayer. This conclusion was 
based on the fact that heat transfer measurements with suspen- 
sions of aluminum powder in argon showed that particle thermal 
properties had no effect on Nusselt number reduction in the 
vertical test section. 

3. Gravity segregation and thermophoresis causing deposition on 
the cool wall of the test section complicated the results in 
the horizontal test section. The Reynolds analogy was not 
considered valid for horizontal flow because the particles did 
not just distort the structure of the turbulent fluid, as they 
did in vertical flow, but actually became an active part of 
the heat transfer mechanism. 

4. The aluminum powder produced unusual deposits. Rock-hard 
deposits appeared on any protuberance into the flow channel. 
The aluminum welded itself to metallic objects by plastically 
deforming on collision with initial non-uniform deposits 
formed by other means. The other means of deposition could 
have been caused by electrostatic attraction, gravity sedi- 
mentation, or thermophoresis. 
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List of Symbols 


A 


heat transfer surface area in equations 2, 3, and 5, ft . 


A. 

1 


inside heat transfer surface area, ft . 

2 

outside heat transfer surface area, ft . 


A 

C 

I 

c 


pg 


pp 


ps 


D 

h, 


ss 


• ^ 

average heat transfer surface area for conduction, ft . 
specific heat capacity of gas, Btulbm 

specific heat capacity of a particle, Btulbm 

specific heat capacity of suspension, Btulbm ^F 

pipe inside diameter, ft- 

-1 -2 -1 

inside heat transfer coefficient, Btuhr ft F 

. . - 1 - 2-1 

outside heat transfer coefficient, Btuhr ft F 
gas thermal conductivity, Btuhr ^ft ^F 

- 1-1 -1 

stainless steel thermal conductivity, Btuhr . ft F 
logarithmic mean temperature difference, F. 


LMTD 

logarit 

-1 

Nu 

h.Dk 

g 

1 g 

-1 

Nu 

h.Dk 

s 

1 g 

Pr 

C Uk 

pg g 


Prandtl number of gas, dimensionless^. 

Q heat transferred, Btuhr 

Re UDv Reynolds number based upon gas properties; 

^ dimensionless. 


thickness of stainless steel tube wall, ft. 


ss 

gas temperature, F. 
inlet suspension temperature, F. 
"^OUT outlet suspension temperature, F. 
T^IN inlet water temperature, F. 
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T 

WOUT 


W 


U 


U 


W 

g 

w 

p 


w 

s 

6 


"n 

p 

V 


outlet water temperature, F. 

average of v inlet and outlet suspension temperatures, F. 

weighted average of wall temperatures in vertical test 
section, F. 

bulk ayerage gas velocity in definition of Reynolds number 
ftsec 

. . - 12-1 

overall heat transfer coeffrcient, Btuhr ft F 
mass flow rate of gas only, Ibmsec ^ or Ibmhr 
mass flow rate of particles only, Ibmsec ^ or Ibmhr 

mass flow rate of suspension, Ibmsec ^ or Ibmhr 

ratio of specific heat of particle to specific heat of gas 
dimensionless. 

W W loading ratio, dimensionless. 

p g 

-1 -1 

gas viscosity, Ibmft sec 

2 -1 

gas kinematic viscosity, ft sec 
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Figure 2 Clean Gas Nusselt Number Versus 
Reynolds Number 
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Beads and 30.0 Micron Aluminum 



Figure 5 Horizontal Nusselt Number Ratio Versus Loading Ratio 
21.6 Micron Glass Beads 
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Figure 6 Horizontal Nusselt Number Ratio Versus Loading Ratio - 36.0 Micron 
Glass Beads and 30.0 Micron Aluminum Powder 
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Figure 7 

Aluminum Deposit on a Thermocouple 
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